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O . The understanding of developed turbulence, a long standing challenge for mathematical physics, 

has entered into the third millennium as an unsolved problem. It poses basic questions con- 
cerning both the behavior of solutions of hydrodynamical equations and the basic principles of 
statistical mechanics of systems out of equilibrium. Although we seem far away from the defi- 
nite answers to these million (or more) dollar questions (Fefferman, 2000), some insight may be 
gained by studying simpler models showing similar behaviors, often themselves not devoided of 
direct physical interest. One of such problems concerns the passive transport of scalar quanti- 
ties such as temperature or tracer 6{t, x) and pollutant or dye density p{t, x) by random flows. 
The flows are described by a simple random ensemble of velocities v(t, x) that is considered 
given. Such approach ignores the back-reaction of the scalar on the velocity dynamics as well as 
VT) • many details of that dynamics. It incorporates, however, the basic phenomenological property 

^ . of velocities exhibiting developed turbulence: their (approximate) statistical scaling 

O 

(N' v(t,x)-v(t,y) ~ |x-y|" (0.1) 

o . 

signaled by the scaling behavior of the expectations of the powers of the velocity differences, 
the so called velocity structure functions, 

> E |v(t,x)-v(t,y)|^ ^ |x-y|^iv (0.2) 
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with (for low A^) not far from the mean fleld theory (Kolmogorov, 1941) value = 
that would correspond to a = 1/3. The time-evolution of the scalar is described by the 
advection-diffusion equation 

dtO + (v ■ V)^ - nV^e = f , (0.3) 

where k is the (molecular) diffusivity constant and /(t,x) denotes a scalar source that one 
may also take random. Given the distributions of the velocities v and of the sources /, 
one inquires about the statistics of solutions 6 of the above equation. As we shall see, very 
simple distributions of v and / lead to stationary statistical states of scalar exhibiting many 
features observed not only in realistic turbulent transport but also in statistics of turbulent 
velocities, in the atmosphere, in aerodynamical tunnels or in sea channels. Those features 
include persistent dissipation of energy, energy cascades from large scales to small ones or 
vice versa and intermittency. Simple models allow a better understanding of the origin of such 
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behaviors and permit to draw some general conclusions, see (Falkovich, Gawedzki & Vergassola, 
2001) for an extensive review and bibliography. 



LECTURE 1 

(transport of scalars by hydrodynamical flow; the role of fluid particle dynamics, single particle diffu- 
sion versus two-particle dispersion) 



1.1. Solution of the advection-diffusion equation, fluid particles 



It is easy to solve the linear equation ( |0.3| ) describing the scalar evolution in d space dimensions 
when the velocity field is sufficiently smooth. 

(i) . For K = and / = 0, the scalar 9{t,yi) is simply constant along the characteristics 

^ = v(.,a.) (1.1) 
as 

which describe the (Lagrangian) trajectories of fluid particles. In other words, the scalar is 
carried by the flow: 

^(t,x) = eis,xt,^is)), (1.2) 

where Xt^^{s) is the Lagrangian trajectory that passes at time t through point x. Note that 
the forward evolution of 6 corresponds to the backward Lagrangian flow. 

(ii) . In the presence of the source /, the scalar is also created or depleted along the trajectory: 

t 

e{t,^) = e{s,xt,^{s)) + J f{a,xt,^{a))da. (1.3) 



(iii). Finally, when k 0, Xt^:x.{s) should be taken as the solution of the stochastic ODE 

dx = \{s,x)ds + V2k, dw (1.4) 

for the Lagrangian trajectories perturbed by the d-dimensional Brownian motion w(s) and the 
right hand side of eq. ( |1.3|) should be averaged over w: 

t 

^(t,x) = E^[e{s,xt^^{s)) + J f{a,xt,^{a))da). (1.5) 



It is then clear that the statistics of 6{t,x.) is determined by the statistics of the (noisy) 
Lagrangian trajectories. The latter may be captured in two steps. First, we may consider the 
transition probabilities for the Markov process given by the solution of the stochastic equation 
( |1.4|) in a fixed velocity field: 

P(v|t, x; s, dy) = E„ 5(y - Xt^^is)) dy . (1.6) 



9. 



Note that the solution (|1.5|) may be rewritten as 



t 

^(t,x) = J P(v|t,x;s,ciy) e{s,y) + J da J P(v|t, x; a, rfy) /(a,y) (1.7) 

s 

or, in the more handy operator notation, as 

9(t) = P(^r\t,s) 9(s) + j* P{w%a) f{a) da . (1.8) 

The latter formuli make sense for k > also in rough (non-Lipschitz) velocities when eq. ( |1.1|) 
does not have unique solutions. Second, in order to take account of the velocity fluctuations, 
we may consider the joint probability distributions of Lagrangian particles 

N 

Vj.iU^^s.dy) = E Y[ P{v\tn,Xn;Sn,dyn) , (1.9) 

n=l 

where t = (t^ . . . ,t^), x = (x^, . . . ,x^), etc. and E stands for the average over the veloc- 
ity ensemble. We shall be interested in the random ensembles of velocities that are station- 
ary, homogeneous and isotropic, i.e. such that the time and space translations and rotations 
v(t, x) v(t + to, i?ox + Xo) are implemented by the measure preserving action of the corre- 
sponding groups on the probability space of velocities. 

1.2. Single-particle diffusion and Richardson dispersion of two particles 

What is the statistical behavior of fluid particles in stationary, homogeneous and isotropic 
turbulent velocities? The rough answers are as follows. For a single particle, one expects a 
diffusive behavior for sufficiently long times. Note that 

t t 
Xo,^it) ~x = J v{a,Xo^^{a))da = Jv^{a,x.)da (1.10) 



for solutions of eq. ( |1.1| ), where v^(cr, x) denotes the the velocity along the Lagrangian trajec- 
tory passing at time zero through x, the so called Lagrangian velocity. For fixed x, {a, x) 
is a stationary process as long as velocities are incompressible. It has zero expectation if the 
similar property holds for the (Eulerian) velocity v((T, x) with which it coincides at time zero. 
If v^(s,x) has temporal correlations that decay fast enough, then the integral in ( |1.10D falls 
under the Central Limit Theorem behaving effectively as a sum of many independent equally 
distributed random variables. As a result, /i~^a;o,x(/^^) tends when fi oo to a Brownian 
motion and 

Pi(0,x;/.t,d(/.-y)) e-'^"^\^,y)dy = j^^^e-^-y^y^'^»'Uy, (1.11) 

i.e. it becomes the transition probability of diffusion. The diffusion constant is given by the 
formula (Taylor, 1921) 

oo 

Do = y'EvjO,x)-vJa,x) rfa. (1.12) 



Still, the decay of temporal correlations of the Lagrangian velocity is a nontrivial fact, see 
(Fannjiang & Papanicolaou, 1996) or (Majda & Kramer, 1999). It does not automatically follow 



from similar property of the (Eulerian) velocity v(a, x) with fixed x and may fail altogether 
in specially prepared velocity ensembles. 



For two particles, the important quantity to study is the 2-particle separation p{t) defined 
as the difference a:;o,x2 (t) ~ ^o,xi (t) of particle positions. It satisfies the equation 

^ = v(s,xo,xi +p(s)) - v(s,a;o,xi(s)) (1.13) 

with the initial condition p(0) = X2 — xi. To get a rough idea about the behavior of the 
particle separation, let us first consider small separations in smooth slowly varying velocities 
where the right hand side may be approximated by linear expression A p{s). This leads to a 
solution 

p{t) - e^*p(0) (1.14) 

with the exponential growth if A has eigenvalues with positive real part (positive Lyapunov 
exponents) that signal the sensitive dependence on initial conditions usually considered as a 
definition of chaos. Notice that although the nearby trajectories separate exponentially in this 
case, the very close trajectories take long time to separate and infinitesimally close ones never 
separate. As the result, the trajectories are still labeled in a continuous way by their initial 
positions. Such fluid particle trajectory behavior pertains to the so called Batchelor regime of 
turbulent flows corresponding to short scales dominated by viscous effects. 



Suppose now that we solve the equation ( |1.13| ) in the regime where the velocity difference 



behaves like p(s)" with a < 1. Such scaling behaviors are observed in the inertial range of 
scales of turbulent flows where viscous and stearring effects are negligible. Passing to the scalar 
version of eq. (|1.13| ), we obtain 



ds ds 



2^-p oc p"+i (1.15) 



and, ignoring again the time and point dependence as well as the statistical fluctuations in the 
proportionality constants, we obtain 

p(t)|i-" ~ P(0)^"" + const, t. (1.16) 

This very rough estimate allows us to expect a power law growth of the 2-particle dispersion 
(the distance between two particles), wiping out the memory of the initial separation. In 
particular, this would mean that inflnitesimally close trajectories would still separate in a flnite 
time, unlike in smooth velocities, leading to a spontaneous randomness in the Lagrangian flow 
at K = 0. Of course, in non-smooth velocities (e.g. in the Holder continuous ones) one should 
not expect existence of the deterministic Lagrangian flow with trajectories labeled by initial 
conditions since the assumptions of the theorem about the uniqueness of solutions of the ODE 
( |1.1D require Lipschitz continuity of velocities in space. That type of situation pertains to 
the behavior of turbulent flows at very high (ideally, inflnite) Reynolds numbers when the 
scaling behavior (|0.1| ) extends down to very small (inflnitesimal) separations crossing over to 
the Lipschitz behavior with a = 1 only at scales where the viscous effects become important 
(of order of fractions of millimeter in the turbulent atmosphere). The super- diffusive behavior 
p^{t) oc corresponding to ( |1.16| ) with the Kolmogorov value a = 1/3 has, indeed, been 



observed phenomenologically for the mean squared 2-particle dispersion. It constitutes the 
content of the flrst quantitative law of developed turbulence formulated in (Richardson, 1926) 



A 



on the bases of experimental data about the separation of meteorological balloons and smoke 
particles. 



The above arguments ignored the temporal dependence of velocity fields which plays an 
important role. Nevertheless, the basic conclusion about the possibility of finite time separation 
of arbitrary close Lagrangian trajectories in spatially rough velocity fields holds even if velocities 
are completely decorrelated (white) in time, as we shall see below. 

The statistics of the separation of two fluid particles may be captured by the relative tran- 
sition amplitudes 



rcl 



ii,Po;s,dp) = y"p2(t,t,xi,xi + Po);s,s,yi,c/(yi + p)) (1.17) 



where the integral is over the final position yi of the first particle. A strong version of the 
Richardson-type super- diffusive behavior for large times may be formulated as the statement 
about the existence of the limit of the rescaled process /i~T^p(/it) or of the limit 



lim Vf{0,Po;fit,d{fiT±^p)). (1.18) 



Note the difference with the expected diffusive scaling ( |1 . 1 1| ) for a single particle. Again 



this type of scaling is not automatically guaranteed by the statistical scaling of the Eulerian 
velocity differences since there are many points where the naive mean-field type arguments may 
go wrong. The super-diffusive behavior of the 2-particle dispersion is, in general, even harder 
to establish than the diffusive behavior of a single particle. 



LECTURE 2 

(Kraichnan ensemble of velocities; 
multi-particle statistics) 



Le Jan-Raimond construction 



of Lagrangian particle processes; 



It is important to have at the disposal a simple model where the ideas about the behavior of 
fluid particles discussed in the first lecture could be tested rigorously. 

2.1. Kraichnan ensemble of velocities 

Such a model has arisen from the work initiated in (Kraichnan, 1968). Kraichnan proposed to 
consider a Gaussian ensemble of velocities decorrelated in time but with the scaling properties 
in space built in. Gaussian ensembles are completely determined by the 1-point and 2-point 
functions. One assumes that the 1-point function of v vanishes and that 

Ev^(t,x)v^(s,y) = 5(t-s)D^^(x-y), (2.1) 

where 

fl-.(x) = / (5"-—) ,, <tk. (2.2) 

lk|<r,-l 

Decomposing 

D'^ = Do 6'^ - rf^^ (x) (2.3) 



with the first term equal to D'^^{0), it is not difficult to see that d^^{x.) — > Do6^^ with 
Do = 0{L^). At short distances 

rf*^(x) = Di[(c/+l)(5^^ |x|2 - 2xV] + 0(|x|^) for |x|<r7, (2.4) 

whereas 

rf^^'(x) — ^ D2[{d-l+^)6'^\x\^ -^x'x^ \^\^-^]. (2.5) 

The limiting scaling formula ( p.5|) approximates well d^^{x.) in the "inertial interval" rj <^ 
x -C L. The scale rj plays the role of the "viscous scale" within which the fractional scaling 
of d^^{x.) is replaced by a quadratic behavior and scale L of the "integral scale" on which 
velocities decorrelate. The matrix d{x.) describes the correlations of the velocity differences. 
E.g. 

E {v\t,x) -v\t,0)){v^{s,^) -v'{s,0)) = 2 6{t-s) d'^{^). (2.6) 

The Kraichnan ensemble of velocities incorporates on the statistical level the scaling properties 
( p.l| ) of velocities in the inertial interval ?7 -C |x — y| <C L. For r] > 0, the Gaussian measure of 
the ensemble is supported by smooth velocities. This is not the case, however, in the limiting 
case 7] = where the viscous scale is set to zero, mimicking the inviscid or inffnite Reynolds 
number ensemble of turbulent velocities. In this case the ensemble measure is supported on 
velocities that in their spatial behavior are Holder continuous with any exponent smaller than 
^/2 whereas the velocities with the Holder exponent bigger than ^/2 have measure zero. 
Of course, in their temporal behavior, the Kraichnan velocities behave as white noise or a 
derivative of the Brownian motion and thus are distributional. We shall have then to modify 
our mean-ffeld arguments adapting it to such a case. 

2.2. Advection-diffusion equation in the Kraichnan velocities 

Let us consider ffrst the finite-dimensional analog of the advection-diffusion equation (|0.3| ), 

e = f3{t)e -aO, (2.7) 

where 6{t) takes values in R-^ and f3{t) is a skew-symmetric matrix (a counterpart of the 
operator — v ■ V) and a a positive symmetric matrix (a counterpart of — kV^). Suppose first 
that /9 is a smooth function of time. The solution of the above linear equation has the form 
9(t) = P{P\t, s) 9{s) with the propagator P{P\t, s) given by the time-ordered exponential 

J (I3(a)-a)da 

P{P\t,s) = T 

oo „ 

= / e-(*-"")"/?(cr„) rfcr„ e-(""-""-i)" ... 

s<ai<a2<...<(T„<t 

... (3{a2)da2 e-("^-"^)'^/5(ai) dai e-^"^-^)" (2.8) 
for t > s. Another way to express the same solution is by a limiting procedure: 

J {0(a)-a) da J (/3W-a) da J {I3{a)-a) da 

P(/3|t,s) = lim e-" e""-i ... e= . (2.9) 

min{(Tm-O-m-l)\0 



Suppose now that is a white-noise Gaussian process with values in antisymmetric matrices 
with mean zero and with the 2-point function 



E/3"^(t)/3^'5(s) = 6{t - s) C"^'^^ . (2.10) 

Now (3 da = dW is a differential of a Brownian marix- valued motion W{t) and the integrals 
in (|2.8|) become stochastic ones. Consequently, we fall into the standard ambiguity with the 
choice of the convention for the latter. On the other hand, the formula ( |2.9| ) still makes sense 
with the convergence taking place in any of the Gaussian process with p < oo, as it is not 
difficult to show. We shall take it as the solution of the equation ( ^.7| ) for the white (3{t). In 
terms of the standard conventions, this corresponds to the Stratonovich prescription in ( p.8|) 
and provides the solution of the Stratonovich stochastic ODE 

de = dWoO - aedt (2.11) 

or, equivalently, of the Ito one: 

de = {dW)e - {a + c)9dt (2.12) 

with c"^ = —1(7"'^''^^ (summation over 7!). The solution of the latter equation, in turn, is 
given by the version of eq. (|2.8| ) with a replaced by (a + c) and the integrals interpreted as 
the Ito stochastic ones and [3 {a) da as dW{a). 



We shall then try to define the the solution of the advection-diffusion equation ( p.3|) for 
white in time velocities f(t,x) dt = (iV(t,x) as given by ( |1.8|) with P{-v\t,s) represented by 

t 

J (-dV(cr)-V+KV2 da) 

P{v\t,s) = T e- 

00 . ^ ^ 

= (-1)" / e dY{an) ■ V e ('^-"'^"-i) 

"■=° S<(7i<(T2<...<cr„<t 

... rfV((T2) ■ Ve("2""i)'^^'rfV((Ti) ■ Ve^"^"')'^^' = ^P„(v|t,s), (2.13) 

n=0 

with the Ito stochastic integrals and k, = k, + i-Do- The operator- valued white noise —dY{t) ■ V 
plays the role of dW{t), the operator kV^ the one of —a, and ^DqV"^ that of — c, as is easy 
to figure out by comparing ( |2.lD and { \i-'S\ ) with ( |2.1CI| ). 

It is not difficult to check that each term of the series is well defined in the action on a 
function, say, from L°°(R'^) and, in its dependence on v, belongs to of the Gaussian 
process with p < 00. The only problem is the convergence of the series that was established 
for K > (k = included!) in (Le Jan & Raimond, 1999). The argument is quite simple. 

N 

First, for / G L°°(R ) and P^j^(v|t, s) = J2 -fn(v|t, s) one shows inductively the a priori 

n=l 

bound: 

E\P^,iv\t,s)f\' < P,it,s)\f\\ (2.14) 

where the O^'^-order term P^it, s) = e But the terms P„(v|t, s) f are orthogonal with 
respect to the scalar product in since all the differentials (iV(crn) are independent in virtue 
of the Ito convention and their expectations vanish. Hence the bound ( |2.14| ) establishes the 
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convergence of the series J2 -fn(v|t, s) f in the norm of the velocity process for any k > 

n 

together with the hmiting bound 

E|P(v|t,s)/r < Poit,s)\f\\ (2.15) 
The continuity of the result in k > may be also established. 



Proof of estimate ( p.l4[ ). For = 0, the bound reduces to the easy estimate \Po{t, s) /p < 



Po(t,s) l/p. Suppose now that ( gUD holds up to N. Note that 



P,^^,(v|t,s)/ = P,{t,s)f - Jp,Mt,(^)dV{a)-VP,{a,s)f. (2.16) 



Squaring and taking expectations, we obtain 

E|P,,^,(v|t,.)/p = \P,{t,s)ff 
t 



+ J rfaE (P,^(v|t,a)®P<^(v|t,(T)) /C2 (Po(a,s)/ ® P,{a,s)f) , (2.17) 

s 

where /C2 = -D*-^(xi — X2)V3,» V^j is the operator acting on functions on R'^ x R'^. Rewriting 
the positive-definite function D^^{x) with the use of the Fourier transform as / D*-^(k) e^^'^ 

. d.-l_ 

with the positive matrix D*-'(k) = J2 -^L(k)A^(k), we may present the second term on the 

0=1 



right hand side of (|2.17| ) as 

t 

[da fdkY^E |P<,(v|t,a)^,k,a|', (2.18) 

s 

where /(j,k,a(x) = e"**^" A^(k) V- [Po(cr, s) /](x). By the inductive hypothesis, this expression 
is bounded by 

t t 

JdaJdkJ2 Poit,^) IUkcI' = Jdajdk P,{t,a)&{k) (v,P„(a,s)7) {VT,{a,s)f) 

t t 
= DoJ daP,{t,a) \VPM,s)ff < 2 k j da P,{t,a) |VP„(a,s)/|' . (2.19) 

s s 

Altogether, we obtain 

t 

E |P,^^,(v|t,s)/|2 < \P,{t,s)jf + 2KjdaPoit,a)\VPM,s)ff 

s 

\ d 

= \Poit,s)fV - J da-{P,Xt,a)\PM,s)ff) = Po(t,.)|/r (2.20) 



ending the inductive proof of eq. ( |2.14| ). □ 



The Chapmann-Kolmogorov chain relation P(v|t, a) P(v|(j, s) = P(v|t, s) for s < a <t 
and the normalization P(v|t, s) 1 = 1 follow easily. Le Jan and Raimond also proved that 
the operators P(v|t, s) preserve positivity. This implies that 

|P(v|t,s)/| < P(v|t,.)|/| < P(v|t,s)||/|L = ll/IL (2.21) 



8 



so that -P(v|t, s) f are (essentially) bounded. We obtain this way a family of Markov transition 
probabilities parametrized by the velocities of the Kraichnan ensemble, hence, also a family of 
Markov processes describing the Lagrangian trajectories (with and without the perturbing 
noise). Taking the N —>■ oo limit in eq. ( p. 161) , we infer that P(v|t, s) satisfies the stochastic 



(anti-)It6 integral equation 

t 

P(v|t,s)/ = P,{t,s)f - Jp{v\t,a)dV{a)-VPM,s)f. (2.22) 

s 

which, upon differentiation over s, gives the stochastic versions of the diffusion-advection 
equation 

ds P(v|t, s)f = P(v|t, s) dV{s) -Vf-K P{v\t, s) VV ds (2.23) 

= P{^r\t,s)odV{s) -Vf - KP{Y\t,s)V'^fds. (2.24) 



2.3. A^-particle processes in Kraichnan velocities 

As follows for example from the relation ( p. 22]) , the expectation of the transition probability 
P{v\t,s) coincides with the first term Pg{t,s) = e^*"^)*^^^ of the series ( |2.13| ). Thus, recalling 
the definition (|1.9|) , we obtain 

V,{t,^;s,dy) = e(*-^)^^'(x,y)rfy. (2.25) 

We infer that in the Kraichnan model, a single fluid particle undergoes diffusion for all times 
with the diffusion constant equal to k, i.e. to the sum of the molecular diffusivity k and 
of the "eddy diffusivity" ^Dq, an effective diffusivity due to the random velocities. Recall 
that Dq = 0{L^) so that the eddy diffusivity is dominated by the integral scale, i.e. by the 
correlation length of the velocities. The virtue of the use of the stochastic Ito integrals in ( p. 131 ) 
was that it made the regularizing role of the eddy diffusion explicit and permitted uniform 
treatment of the cases with k > and with k = 0. 

In order to study the statistics of particles in the Kraichnan model, we have to analyze 
the joint transition probabilities V^{t,x,;s,dy), see ( p..9|) . If fact, it is enough to look at 
their equal-time versions Vj^{t,x^; s,dy). The stochastic (anti-)It6 equation ( p. 231) implies the 
relation 

^ [ P,(t,x;.,rfy)/(y) = - E / V^t^TL; s,dy) /^^^(yn - y„.) V,. V^^/(y) 

-Y. fK{t,^;s,dy)RVlf{y) (2.26) 

from which one deduces that 

V^{t,2L;s,dy) = e(*-^)^-(x,y) rfy (2.27) 
for the second order differential operator 

= ^E^^^(x.-xjV.,V,.^ + '^EVL- (2.28) 



It follows that in the Kraichnan model the transition probabilities Vj^{t,x,', s,dy) still form a 
Markov family (this is due to the temporal decorrelation of velocities). The probabilities of 
relative separations of particles 

Pl'\t,^;s,dy) = y'p^(t,x;s,rf(yi + y),...,ci(y^ + y), (2.29) 

with the integral over the translations y, are given by the exponential function of the operators 
that are restrictions of Aij^^s to the translation-invariant sector: 

^jv = - E {d'\^n - x^) + 2/s:) V,. V^.^ (2.30) 

n<m 

with (i*-' (x) given by (|2.3| ). As we see, in the Kraichnan velocities, fluid particles undergo in 
their relative motion an effective diffusion process with the configuration-dependent diffusivity. 

In particular for the probability distribution ( |1.17] ) of the relative separation of two fluid 
particles, we obtain 

Vf{t,p^-sAp) = e(*-^)*^Hpo,P)rfP, (2.31) 

where 

M2 = S^{p)V,Vj + 2/tV^ (2.32) 

This operator commutes with the action of the rotation group in L^(R'^) and reduces in 
the action on functions carrying irreducible representations of SO{d) (labeled by the angular 
momentum i = 0,1, . . .) to a second order differential operator in the radial variable p = \p\. 
In particular, the probability distribution V^it, po; s,dp) of the 2-particle dispersion is given 
by the exponential function of the restriction M|^° of operator M2 to the rotation-invariant 
sector with i = 0. 

The above relations will allow us to analyze in the next two lectures the properties of the 
fluid particles in the Kraichnan ensemble of velocities. 



LECTURE 3 

(fluid particles and advection of scalar in the Batchelor regime of the Kraichnan model: random chaos) 

3.1. Separation of close particles in smooth Kraichnan velocities 

In the region where the distances between the fluid particles are much smaller then the viscous 
scale Tj, i.e. in the Batchelor regime, we may approximate d^^ as in ( p.4|) . Upon dropping 
the 4*'^-order terms (which, strictly speaking pertains to the behavior of infinitesimally close 
trajectories) and upon setting k = 0, the generators of the A^-particle processes become 

= -Di E ((^n - ^-)'Vx„ ■ + 2 (xl, - xU{x\ - x^„)V.; V,.^) . (3.1) 
In particular, for N = 2, 

M2 = D,[p^V^ - 2pyV,V,] (3.2) 



in 



in terms of the separation variable p. In the rotational invariant sector, 

= D^{d~l)p-'^+^dpp'^+^dp. (3.3) 

From the latter expression, one infers easily the probability distribution of the 2-particle dis- 
persion: 



VMPo:t,dp) = -^^=^^^ exp[- j^5^(^(ln(^)-Di(rf-l)rft) ] ^. (3.4) 



^ p^r. — ~ nJd-^^dt' 

^4iTDi{d-l)t 

As we see, the logarithm of the 2-particle dispersion grows linearly with the rate (the Lyapunov 
exponent) X = Di{d — l)d > 0. The system exhibits the exponential separation of trajectories, 
i.e. is chaotic. It is easy to see that 

lim V^{0,po;t,dp) = 5{p)dp, (3.5) 

which is a consequence of the existence of deterministic Lagrangian trajectories in smooth 
velocities: it signals that the Markov process with the transition probabilities P(v|t, x; s, (iy) in 
a fixed velocity v realization concentrates (for k = 0) on deterministic trajectories determined 
by the initial conditions. 

It is not difficult to understand the origin of the positivity of the Lyapunov exponent in the 
smooth Kraichnan velocities. For such velocities, the equation ( p,.13[ ) for very close trajectory 
separation may be approximated by 



dp 

ds 



p- Vv(s, xo,xi(s)). (3.6) 



In the Kraichnan model, 

EVy{t,^)Vy{s,y) = 6{t-s)V,V,d'^{^-y) (3.7) 

which is independent of (x — y) in the quadratic approximation to (i*''(x). In other words, 
in this approximation, Vj.t'*(t,x) = '-j^^{t) where 7*'^(t) is a (Gaussian) white noise with the 
values in the real traceless matrices, with mean zero and covariance 



E f\t) -f^\s) = 2Di 6{t - s) [{d + 1) 6'' 5^^ - 5'^ 5'^ - 5'^ 5'^) . (3.8) 

Eq. (|3.6|) becomes the stochastic differential equation 

dp = dVop (3.9) 

where r(s) is a Brownian motion on traceless matrices with dr{s) = 7(5) ds. The solution is 
given by 

p{s) = Gt{s)p{t), (3.10) 

where Gt{s) is a Brownian motion on the group SL{d) of real unimodular matrices satisfying 
Gt{t) = 1. It follows that in the quadratic approximation to d^^{'x), 

J P(v|t, xi; s, dy,) P(v|t, X2; s, dy^) /(ys - Yi) = / {Gt{s)p,) , (3.11) 



n 



where po = ^2 — xi. Taking the averages of the last identity, we infer that 

(pf(t,.)/)(Po) = Bf{G,{s)p,) . (3.12) 

It is not difficult to find the generator of the Brownian motion Gt{s). Consider the natural 
actions of SL{d) and of its subgroup SO{d) in L^(R'^). Their infinitesimal generators are 

= (-p^v.+p^V.) , (3.13) 
respectively. A simple algebra shows that 

M2 = Di[dH^ ~ {d+l)J^] (3.14) 
where and are the quadratic Casimirs of SL{d) and of SO{d) : 

= H'^ W\ J2 = -i(/^y. (3.15) 

The right hand side of ( p.l4| ) (with the Casimirs interpreted as those of the left regular action) 
gives the generator of the Brownian motion Gt{s) on SL{d). It is well known that such 
Brownian motions, that may be thought of as continuous products of random, independent, 
identically distributed matrices in SL{d), lead to positive Lyapunov exponents (a continuous 
version of the Furstenberg-Kesten or Oseledets Theorems, see e.g. (Arnold, 1998)). 

For the A^-point operator given by ( |3.1| ), we similarly obtain: 

= D,[dHl ^ - id + l)JlJ (3.16) 

where _ and are the quadratic Casimirs of the diagonal actions of SL{d) and of 

SO{d) on the (A^ — 1) separation variables, e.g. on p„ = x„ — xi. An alternative expression 
for is 

= D,[dGl^^ + ^^A{A+{N-l)d) - {d+l)JlJ (3.17) 

where G^ .^ is the quadratic Casimir of the action of SL{N — 1) on the index n of pl^ with 
the generators G""* = — p^ Vpj^ + j^^y ^"'"^ S ^p^xj, and A = ^p^x^ is the generator of 

the overall dilations. Consider a function / of A^ — 1 separations that depends only on the 
volume p = ^det„„(p„ ■ p„) spanned by A^ — 1 separation vectors (for A^ — 1 < d). Since 
p is SL{N — 1) and S'0((i)-invariant, it follows from ( p.l7| ) that M^f is still a function of p 
only and that 

(MJ) (p) = D, ^ A(A + {N- l)d) f{p) 

= D,{N-l)id-N+l) (p-'^+i Spp'^+i dp /(p) . (3.18) 

Comparing to ( |3.3| ) and ( p.4| ), we infer that the logarithm of p grows linearly with the rate 
Di{N — l){d — N + l)d. By definition, the latter gives the sum Ai + . . . + X^_-^ of the (A^ — 1) 
biggest Lyapunov exponents so that 

\n = Di{d-2N+l)d. (3.19) 
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In the smooth d-dimensional Kraichnan flow the d Lyapunov exponents are equidistant, with 
the sum equal to zero (as required by incompressibility) . The relations of the multi-trajectory 
processes to the harmonic analysis on the groups SL{N) were observed in (Shraiman & Siggia, 
1995 and 1996), 



3.2 Scalar statistics in the Batchelor regime 



One may easily control the stationary state of the scalar developing under the forced transport 
of Kraichnan velocities in the Batchelor regime. For white in time velocities, in agreement with 
the preceding discussion, the solution of the forced transport problem is given by eq. ( p..8| ) with 
the evolution operators P(v|t, s) as constructed above. For simplicity, let us take the zero 
initial data for the scalar: 6{s) = 0. Then the characteristic function of the time t probability 
distribution of the scalar, 

t 

ijdcrjd^ /i(x) P(v|t,x;cr,dy) /(cr,y) 

$(/i|t,s) = ^; eM = Ee ^ . (3.20) 

Let us take the source /(s,y) to be also a Gaussian process, independent of velocities, with 
mean zero and covariance 

E fit, x) /(s, y) = 5it - s) x(x - y) , (3.21) 

where x is a smooth, positive-definite, decaying function on R''. For such a forcing, the 
expectation over the source / of the scalar in ( p. 20]) may be easily calculated, leaving us with 
the expectation over the velocity ensemble: 

t 

-\ j da j dxi dx2 /i(xi)h.(x2) P{v|t,xi; (T,dyi) P(v|t,X2; CT,c(y2) x(y2-yi) 

^h\t,s) = Ee s . (3.22) 

In the Batchelor regime, we may use the relation ( |3.11| ) to rewrite the latter expectation in 
terms of the one over the Brownian motion Gt{(y) on SL{d): 

t t 

-i/dcr/ dxidxa Mxi)/i(x2) x(Gt((T)(xi-X2)) -/ Vh(Gt(a)) da 

$(/i|t,s) = E e ^ = E e s ^ (3.23) 

where VhiG) = h 1 1 d'Xi dii.2 ^(xi) /i(x2) x (G'(xi — X2)) is a positive potential on group SL{d). 
With the use of the Feynman-Kac formula, eq. ( |3.23| ) may be expressed as 

Hh\t,s)= I e(*-^)['^^'-('^+^)-^'-^''](l,G') rfC 

SL{d) 

= 1- r^da (e'^['^^'-("+i)^'-^'']l^,)(l), (3.24) 



where the last equality follows by integration by parts. It is then easy to show that the limit 
$(/i) = lim ^(h\t,s) exists and is a characteristic function of a probability measure supported 

s— >— 00 

e.g. on the space iS'(R'^) of tempered distributions. The latter describes the stationary state of 
the forced scalar which, indeed, is not supported by smooth scalar configurations, as reflected 
by logarithmic singularities at coinciding points of the scalar A^-point functions. For x a^id h 
rotation-invariant, potential Vh is left and right S'0((i)-invariant. On SO{d)\SL{d)/ SO{d), 
the operator [dH"^ — {d + 1) J^] reduces the the Calogero-Sutherland integrable Hamiltonian 



and the analysis of the distribution of the random variable / /i(x) ^(x) = 6{h) reduces to the 
analysis of the property of its perturbation by the potential proportional to Vh- For example, 
the exponential rate of decay of the probability density function of 6{h) is related to the bound 
state energy of such a perturbation with a negative coefficient. The details may be found in 
(Bernard, Gawedzki & Kupiainen, 1998). 



LECTURE 4 

(fluid particles in non-smooth Kraichnan velocity fields; breakdown of deterministic Lagrangian flow; 
intrinsic stochasticity versus particle aggregation) 



4.1. Separation of close particles in non-smooth Kraichnan velocities 

In the limiting case rj oi the Kraichnan model, the Lagrangian trajectories exhibit even 
more dramatic behavior. For rj = and L = 0, when d^^{x.) takes the scaling form of eq. ( [^.5| ), 
and for k = 0, 

Mf° = D^id-l) p-^+^dpp'^-^^^dp. (4.1) 

The probability distribution of the 2-particle dispersion 

VMPo;t,dp) = e'''t\p^^p)dp (4.2) 

may be studied for example using the spectral decomposition of Ml^""'^ . The long-time- large- 
distance asymptotics of the dispersion follows from the rescaling property 

P,(o,po;/it,rf(//^'"«V)) = n(o,/i-i/(2-«)po;t,rfp) (4.3) 

and the easily shown relation 

lim P2(0,po;t,f^p) oc p'^-ir''/(2-0 e— t-p'-V* rfp. (4.4) 

po^O 

In particular, one obtains the Richardson dispersion law in the form 

Ep(t)2 = I p'vMpo;t,dp) = Oit"^') 

for large times, with an exact proportionality for all times in the limit po — ^ 0. Such a behavior 
may be predicted by solving the modified version 

c/p2 oc rfwp«/^ + ^ (4.5) 

of the naive mean field type equation (|1.15|) , where w{t) is the Brownian motion introduced 
to account for the temporal decorrelation of Kraichnan velocities. In particular, the growth of 
p^ proportional to is obtained for ^ = 4/3. 

Note that the limit (|4.4| ) is absolutely continuous with respect to the Lebesque measure dp, 
in contrast to what happens in the Batchelor regime, see eq. ( p.5|) . Such a property excludes the 
concentration of the transition probabilities P(v|t, x; s, dy) at a single (v-dependent) point y 
and signifies the breakdown of the deterministic Lagrangian fiow with fiuid particle trajectories 
determined by their initial positions in fixed velocity fields. Instead, the Lagrangian trajectories 



U 



form a genuinely stochastic process already in a fixed velocity realization, as predicted above 
and illustrated in the figure below: 






space 



That Markov process, may be obtained by first adding the regularizing noise to the trajec- 
tory equation, as in (|1.4|) , and then turning it off, or also directly, by the Le Jan-Raimond 
construction of the n = transition probabilities P(v|t, x; s, dy). 

4.2. The role of compressibility 

Up to now, we have considered the Kraichnan ensemble with incompressible velocities, the 
property guaranteed by the presence of the transverse projector (S^^ — ^p^) in the Fourier rep- 
resentation ( |2.2| ) of the velocity 2-point function. In order to study the effects of compressibility, 
one may introduce another parameter besides the roughness exponent ^, the compressibility de- 
gree p. We shall do it by replacing the transverse projector in (|2.2|) by [{l — p)S^^ + {pd~l)^^]. 
The value p = corresponds to the incompressible case whereas for p = 1 almost all veloci- 
ties are gradients, with the intermediate values of p interpolating between the two cases. The 
preceding constructions, in particular the Le Jan-Raimond one, carry over to the case with 
non-zero p. The scaling form of the generator of the 2-particle dispersion process becomes 
now 

M^=^ = D2{d-1) p^-^dpp^dp (4.6) 
with a = ji^^ — 1- The definition of the semigroup e*^2 ° requires a choice of the boundary 



condition for the operator Mg^" at p = 0. Such a choice is automatically assured by con- 
sidering first the k > case and then sending k to zero. This limiting procedure selects for 
p < d/^"^, i.e. for weak compressibility, the eigenfunctions of M^^" behaving like 0{1) at p = 
whereas for strong compressibility p > d/^'^ the eigenfunctions that behave like 0{p^^°') are 
chosen. The two choices result in very different probability distributions V^iO, po;t,dp) of the 
2-particle dispersion and, in consequence, to dramatically different Lagrangian flows, as first 



noticed in (Gawedzki & Vergassola, 2000). For p < d/^ , in a simple generalization of ( [4.4|) , 



lim V2{0,po;t,dp) oc p"-€t(€-i-")/(2-0 e-^^^^t-^'^Vt (4.7) 

from which the conclusions about the the intrinsic stochasticity of the Lagrangian flow may be 
drawn the same way as for the incompressible case. For p > ^/d"^, however, 

V,{0,po;t,dp) = V''^{0,po;t,dp) + p{t,po)5{p)dp (4.8) 

with V'^'^^{0, po;t,dp) absolutely continuous with respect to the Lebesque measure dp and 
p(t,po) > 0. When po — the regular part tends to zero and p(t,po) tends to 1 so that one 
recovers the behavior ( p.5|) . The latter indicates that the Lagrangian flow is deterministic, with 
trajectories determined by the initial condition in fixed velocity realizations. The presence of 
the term proportional to the delta- function for po > signals, however, that, with positive 
i-dependent probability, the trajectories starting at different initial points collapse together by 
time t, as illustrated in the figure on the next page. Such a non-conventional behavior of 
trajectories is again possible since the theorem about the existence and unicity of trajectories 
fails for non-Lipschitz velocities. The competition between the tendency of the trajectories in 
such velocities to separate explosively and their trapping in the regions of strong compression 
is won by the first trend for weak compressibility and by the second one for the strong one. 

4.3. Persistence of scalar dissipation 

The unusual behaviors of Lagrangian trajectories in rough Kraichnan velocities are the source 
of phenomena crucial for the scalar transport. One of the implications of the advection-diffusion 
equation ( p.3| ) in the incompressible velocities is the scalar "energy" / 6^ balance 

^ J e{t, x)2 dx = -2k J {V0{t, x))2 dx + 2 y e{t, x) f{t, x) dx (4.9) 

with the first term on the right hand describing the dissipation and the second one the injection 
of scalar energy by the sources. One may naively expect that in the absence of sources the 
scalar energy is conserved in the limit k ^ 0. Let us examine this question more closely. 
Recall that the evolution of scalar is related by eq. ( |1.7| ) to the transition probabilities of the 
(noisy) Lagrangian trajectories. For vanishing sources, one obtains the identity 

1 rfx 1 P(v|t, x; s, dy) [e{s, y) - e{t, x)]^ = J e{s, dy - J 9{t, x)^ rfx (4.10) 



which follows by developing the square and using eq. ( p..7|) , the normalization of the measure 
P(v|t, x; s, (iy) and the symmetry of the operator P(v|t,s). The left hand side is obviously 
non-negative which implies that the scalar energy cannot grow. It vanishes if and only if 
9{s,y) = 0{t,x) on the support of P(v|t, x;s,y) for each x. This happens for arbitrary 
9{s, y) if and only if P(v|)f:, x; s, dy) is supported by exactly one (x-dependent) point, i.e. if the 
Lagrangian trajectories are determined by their single-time positions. The intrinsic stochasticity 



Ifi 



space 



of the Lagrangian flow results then in the persistence of scalar energy dissipation in the limit 

The same effect may be seen in averaged quantities, also in compressible Kraichnan ve- 
locities. Suppose that at the initial time s we are given an homogeneous, isotropic scalar 
distribution with the 2-point function 



Then, if the initial scalar distribution is independent of velocities and in the absence of sources, 
the scalar 2-point function at the later time t is given by 



Ee{s,y)9{s,y + p) 



(4.11) 



F,{t,po) = B9{t,x)9{t,^ + Po) = I F,{s,p)V,{t,po;s,dp) 



(4.12) 



In particular, the mean scalar energy density 




(4.13) 
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Since F^{s,p) < F^{s,0) by the Schwartz inequahty, it follows that the mean scalar energy 
density is non-increasing. It is conserved in general only if lim VoU, po; s, dp) = 6(p) dp), i.e. if 

Po-*0 

the trajectories are uniquely determined by their final positions. This is the case for strongly 
compressible Kraichnan velocities with p > d/^"^. 



LECTURE 5 

(forced scalar advection; intermittent direct cascade with persistent dissipation and zero mode domi- 
nance versus non-intermittent inverse cascade) 

When scalar is forced with the random Gaussian source characterized by the isotropic 2-point 
function (|3.21|) and if the forcing is independent of velocities and of the initial distribution of 



the scalar, then the scalar 2-point function evolves according to the relation 

t 

F^{t,po) = J F^{s,p) V^{t,pQ]s,dp) + j da j xip) V^i^, po; a, dp) , (5.1) 



which solves the differential equation 



|F,(t,p) = Mt"F,it,p) + x(p). (5.2) 

Setting p = in the latter identity, we obtain the averaged scalar energy balance with 
M|"° ^^(t, p)|p^p representing the mean dissipation rate e and x(0) the mean injection rate 
of the scalar energy per unit volume. 

The first term in the solution (|5.1|) decays with time if the initial 2-point function F^{s,p) 
decays in space. The second term may be interpreted as the lapse of time between moments 
s and t when the two trajectories starting at distance po stay in the region of sizable values 
of the source covariance x- Its asymptotic behavior depends crucially on what two Lagrangian 
trajectories do at long times. 

5.1. Direct versus inverse scalar cascades 

In the weakly compressible regime p < d/^"^, the trajectories separate to a fixed distance in 
a finite time. Nevertheless, the average time spent by two trajectories within the range of x 
is finite only for p < ^—^^ (or a < 1). In that case, the scalar 2-point function reaches the 
stationary form independent of its initial value 

t oo pi 

F^{p) = j da j xip) V,{t,po;a,dp) = J Pi" dpi J pt^ x{p2) dp2 , (5.3) 

-oo p 

where the last equality holds for vanishing rj, L and k. If ^-^^ l£ p < d/^'^ then, although 
the trajectories separate to a fixed distance in a finite time, with positive probability they 
revisit smaller separations. As a result, the average time they spend within the range of x 
and, consequently, the scalar 2-point function diverge when t — oo. What still reaches the 
stationary form, however, is the 2-point scalar structure function 

p pi 

E[^(t,x)-0(t,x + p)]2 ^ 5,(p) ^ -1— J p-- dpi J pt^ X{P2) dp2 (5.4) 
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which in the hmit exhibits the scahng behavior ~ ^ for small p. In the stationary state, 
the scalar energy balance reduces to the identity 

Mt'FM\..o + X(0) = (5.5) 

expressing the fact that the dissipation and injection balance each other. In particular, the 
mean dissipation rate e, equal for k > to 2k E (V6')^, is /t-independent and does not vanish 
when K — >^ 0, the phenomenon called the dissipative anomaly. The anomaly, accompanied 
by the direct scalar energy cascade from the scale on which x{p) decays (where it is injected) 
to smaller and smaller scales, is another manifestation of the persistence of the scalar energy 
dissipation. It is assured by the explosive separation of the Lagrangian trajectories in the whole 
p < d/^"^ range. 

In the strongly compressible regime p > d/^'^, the scalar 2-point function does not stabilize 
but has a constant contribution growing linearly in time with the rate equal to x(0) the 
K — *• limit. No dissipation persists in this limit due to the deterministic character of the 
Lagrangian trajectories. The injected scalar energy ultimately condenses in the constant mode 
in the process of the inverse cascade towards longer and longer distances. The 2-point 
structure function of the scalar, however, reaches the stationary form 

p oo 

= I^b) /p^^Pi / PrHxiO) - xiP2)) dp2 . (5.6) 

pi 

with the scaling behaviors oc p^^" for small p and oc p^^'» for large p. 
5.2. Zero mode scenario of intermittency 

In the presence of stationary, Gaussian, time decorrelated sources, the higher-point equal- 

N 

time correlation functions E H G{t,^n) = -^jv('''!^) of scalar satisfy the evolution equations 

n=l 

generalizing ( |5.2| ): 

^F^(t,x) = M^F^(t,x) + (F^_,® x)(t,x), (5.7) 

where 

iFN-2® X){t,^„- ■ ■ = ^ F^_2(t,x,,.^.^.,x^) X(xn -xj . (5.8) 

n<m n m 

The above relations are solved inductively by the expressions 

t 

F^(t,x) = J F^{s,y)V^{t,^;s,dy) + J da J {F^_, ® x){c^,y) V^{t,^;a,dy) , (5.9) 

s 

compare to ( p . 1| ) . One expects that for sufficiently weak compressibility those solutions reach 
stationary form (x) vanishing for odd and given for even by the second term on the 
right hand side with s = — oo. This has been established rigorously for p = and all k, > 
in (Hakulinen, 2002). 

An important question concerns the behavior at small p of the k = stationary state 
scalar structure functions 

s,{t,p) = E [e{^)-e{^ + p)f (5.10) 
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with even N. They are given by special combinations of the correlation functions (x) . Naive 
dimensional predictions based on eq. ( ^.71 ) would suggested the behavior 5*2^^ (p) oc p^(2-5) gj^ce 
operators M^v have dimension length^~^. This agrees with the scaling of the 2-point function 
described above and would automatically hold for the higher structure functions if the scalar 
differences [^(t, x) — ^(t, x + p)] were normally distributed. The dimensional predictions, first 
postulated in (Obukhov, 1949) and (Corrsin, 1951), are in contradiction with experimental 
turbulent advection measurements which indicate scaling of scalar structure functions with 
exponents that grow slower than linearly with A^. Such behavior of the scaling exponents 
signals more frequent appearance than in normal distributions of large fluctuations of the scalar 
differences at small separations, the phenomenon called scalar intermittency. 

It had been suggested in (Kraichnan, 1994) that the higher scalar structure functions of the 
Kraichnan model exhibit non-dimensional scaling. It was subsequently realized in (Shraiman & 
Siggia, 1995), (Gaw§dzki & Kupiainen, 1995) and (Chertkov, Falkovich, Kolokolov & Lebedev, 
1995) that, in the incompressible model, the k = higher point functions Fj^{2^) are dominated 
at short distances by the contributions from the scaling zero modes <^^{x,) of the operators 
satisfying 

M^¥.^(x) = 0, ^^(Ax) = A^-^^(x). (5.11) 

The scaling dimensions Cjv of such modes are not constraint by the dimensional analysis but 
are accessible perturbatively or numerically (note that the modes annihilated by drop out 



in the stationary version of eq. (5/7)). The perturbative calculation of such modes gives for 
even N and general p 

In particular, the zero modes dominate at short distances the structure functions so that 

SM p^-. (5.13) 

Such non-dimensional scaling signals the short distance intermittency of scalar advected by the 
weakly compressible Kraichnan model velocities. 

The zero mode dominance of the stationary scalar higher-point functions has been ex- 
hibited by the perturbative analysis of the Green functions of operators M^^ around ^ = 0, 
d = oo and ^ = 2. The numerical results gave a similar picture for all values of ^, see the 
figure on the next page representing the values of the 4-point function anomalous exponent 
(2^2 ~ C4) obtained by Frisch, Mazzino Noullez & Vergassola (1999) in numerical simulations 
of the three-dimensional (circles) and 2-dimensional (stars) incompressible Kraichnan model. 

What is the physical meaning of the zero modes of the operators that dominate the 
short- distance asymptotics of the scalar A^-point functions? They are statistically conserved 
modes of the effective diffusion of Lagrangian trajectories with generators M^. Indeed, the 
mean value of a translationally-invariant scaling function ip^ (x) of scaling dimension a, viewed 
as a function of time t positions of Lagrangian trajectories, is 

/V^jv(y)n(0,x;t,rfy) (5.14) 

which, for generic ipN^ grows dimensionally as 0{t'^^^'^'^^) for large t refiecting the super- 
diffusive growth oc t^/(2~?) of the distances between the trajectories. But if -j/^^ = ^ 
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zero mode of M^^ then the above expectation is conserved in time (such conserved modes are 
accompanied by descendent ones whose Lagrangian averages grow slower than dimensionally, 
see (Bernard, Gaw§dzki & Kupiainen, 1997). 
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5.3. Non-intermittency of the inverse cascade 

In the strongly compressible phase with the inverse cascade of scalar energy, the behavior of 
the higher structure functions is different. In fact, only the lower ones stabilize, but the ones 
that do, scale normally on large distances. In this regime one can find exactly the stationary 
form of the probability density function of the scalar difference: 

E 5{& - p-M/2 [^(x) - ^(x + p)]) oc [x(0) + const. (5.15) 

at large distances. Its scaling form indicates that there is no intermittcncy in the inverse cascade 
of the scalar (the deviation from the normal distribution is scale- independent). For small p, 
however, all the stabihzing structure functions scale as signaling an extreme short distance 
intermittency. 
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CONCLUSIONS 



As we have seen, the transport of a scalar quantity by velocities distributed according to the 
Kraichnan ensemble shows two different phases characterized by different direction of the scalar 
energy cascades and different degrees of intermittency. The phase transition occurs at the value 
p = of the compressibility degree, where the behavior of the Lagrangian trajectories changes 
drastically from the explosive separation to the implosive aggregation. These two phases are 
somewhat reminiscent of the behavior of the three-dimensional versus two-dimensional devel- 
oped turbulence. That suggests that one should put more stress on the Lagrangian methods 
in studying the latter, not quite a new lesson, see e.g. (Pope, 1994), by with the new twist 
pointing to the importance of the intrinsically stochastic character of the Lagrangian flow at 
extreme Reynolds numbers. Of course, the Navier-Stokes and the Euler equations, unlike the 
scalar advection one, are non-linear, a difference that, certainly, is far from being minor. Also, 
they describe velocity fields that are temporally correlated and transformed when carried along 
their own Lagrangian trajectories. Besides, due to pressure, there are non-local interactions 
present. Some of those effects, however, may be studied already in synthetic velocity ensembles. 
It seems that the investigation of such ensembles has a potential to teach us important lessons 
that have to be mastered on the way to an understanding of fully developed turbulence. 
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